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ABSTRACT 


Transverse  projectile  impacts  on  the  ends  of  free-free  beams  of 
circular  cross  section  produced  evidence  of  hinges  in  the  beams  and  strain 
distributions  which  could  be  used  for  comparisons  with  a  rigid-plastic 
approach  to  the  dynamic  response  of  structures  and  for  comparisons  with  a 
finite  element  computer  code  solution.  The  rigid-plastic  approach 
describes  the  mechanism  of  bending  and  hinge  development  and  motion  correctly 
but  the  simplifications  inherent  in  the  model  prevent  it  freon  accurately 
predicting  the  detail.  The  computer  code  EPIC-2  gives  extremely  good 
simulations  of  the  problem. 
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a,  aQ,  a1 
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linear  acceleration 

D 
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beam  diameter 

1 

= 

beam  length 
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mass  per  unit  length 

Mo 
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limit  moment  of  beam 

Mx 
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bending  moment  at  position  x 

p 

= 

load 

Vx 
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shear  force  at  position  x 

X,  Xh 
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position,  hinge  position 

Y 
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yield  stress  of  beam 

a,  V  a1 
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angular  acceleration 

(Greek,  9 
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A 

3 

Mathematical  Symbol  ~  small  change 
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e 

3 

Strain 

(Greek,  L*C, , 
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BENDING  OF  BEAMS  SUBJECTED  TO  TRANSVERSE 

IMPACTS 


1 .  INTRODUCTION 


When  a  slender  beam  is  subjected  to  a  transverse  impulsive  force 
bending  moments  result  due  to  the  inertia  of  elements  of  the  beam  away  from 
the  site  of  loading.  If  the  pulse  is  of  sufficient  magnitude  then  yielding  or 
fracture  may  occur  at  the  position  of  maximum  bending  moment.  Zener  and 
Peterson  [1]  demonstrated  that  when  an  impact  occurs  at  the  end  of  a  freely 
suspended  beam  the  point  of  maximum  bending  moment  occurs  at  a  distance  one 
third  of  the  beam  length  from  the  impact  site.  Hie  problem  of  impact  at  the 
centre  of  free-free  beams  of  either  infinite  [2]  or  finite  [3]  dimensions  has 
been  given  more  attention;  under  these  conditions  yielding  occurs  initially  at 
the  impact  site  about  which  two  halves  of  the  beam  rotate.  Symonds  [4]  and 
Symonds  and  Leth  (5]  examined  how  the  detail  of  the  loading  and  unloading 
cycle  affects  the  nature  and  position  of  plastic  hinges  developed  in  centrally 
struck  beams,  and  the  influences  of  transverse  shear  and  rotary  inertia 
effects  have  been  considered  by  Karunes  and  Onat  [6]  Symonds  [7]  and  Jones  and 
Gomes  de  Oliveira  [8].  The  main  aspects  of  the  theory  have  been  reviewed  in 
Johnson  [9]  and  Goldsmith  [10]. 


The  present  work  concentrates  on  transverse  impacts  on  one  end  of 
freely  suspended  beams,  as  it  was  found  that  under  these  conditions  plastic 
hinges  can  be  more  easily  produced  away  from  the  impact  site  than  with  central 
impacts.  Ihe  work  has  direct  relevance  to  the  behaviour  of  slender  rod  shaped 
projectiles  which  break-up  [1]  or  bend  [11]  at  specific  points  along  the 
length  of  the  rod  during  impact.  The  results  of  experiments  are  compared  with 
calculations  based  on  the  rigid-plastic  theories  described  above,  and  with 
simulations  using  the  finite  element  computer  code  EPIC-2  [12,13]. 


2.  EXPERIMENT 

Beams  of  circular  section  were  lightly  glued  to  a  support  at  one  end 
and  struck  transversely  with  a  projectile  on  the  other  end.  The  bond  to  the 
support  was  very  weak  so  that  the  beams  readily  separated  and  were  caught 
approximately  0.5  m  away  in  a  soft  medium.  The  mild  steel  impactors  were 
4.76  mm  diameter  cylinders  of  mass  1.65  g,  and  on  impact  at  775  m  s'  they 
mushroomed,  causing  some  local  deformation  to  the  end  of  the  beam.  The  impact 
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conditions  were  estimated  to  be  sufficient  to  exceed  the  shear  yield  stress  on 
the  beams.  All  post-impact  measurements  on  the  beams  were  taken  from  the 
centre  of  the  impact  site  which  was  considered  as  the  end  of  the  beam  for  both 
the  actual  value  of  length  to  diameter  ratio  quoted  and  for  the  strain 
distributions. 

The  beams  were  circular  in  section,  of  diameter  6. 35  mm,  and  of 
nominal  length  to  diameter  ratios  8,  12,  16  and  24.  They  were  manufactured 
from  a  heat  treatable  steel  (Vibrac  R4)  and  used  at  various  hardness  levels 
from  355  HV10  to  575  HV10. 


After  impact  the  beam  profile  was  projected  onto  a  screen  and 
traced,  and  the  drawing  divided  into  a  series  of  segments  along  the  length. 
Using  the  end  of  the  beam  away  from  impact  to  set  a  zero,  the  angle  of  the 
tangent  to  the  beam  was  measured  as  a  function  of  position  along  the  beam  on 
both  the  tension  and  compression  sides.  Assuming  that  the  deformation  is 
purely  bending,  the  strain  s  in  the  outer  fibres  of  a  segment  of  beam  of 
width  A1  is  given  in  terms  of  the  diameter  of  the  beam,  D,  and  the  change  in 
tangent  angle,  A6  ,  across  the  segment  by 


DA9 

2A1 


(1  ) 


The  strain  distribution  as  a  function  of  position  along  the  beam  can  thus  be 
calculated.  The  results  from  compression  and  tension  sides  of  the  beam  were 
similar,  as  expected  for  bending  deformation,  and  hence  only  results  from  the 
tension  side  are  reported. 


3.  RESULTS 

Figure  1  shows  the  profiles  of  impacted  beams  of  two  strength  levels 
and  of  various  length  to  diameter  ratios.  In  each  beam  except  the  shortest, 
in  which  it  would  be  difficult  to  distinguish  two  bends  because  of  their 
likely  proximity,  there  are. two  major  bends;  one  is  approximately  the  same 
distance  from  the  impact  site  in  each  beam  regardless  of  strength  or  length  to 
diameter  ratio  and  the  other  is  approximately  half  way  along  each  beam.  Fig. 

2  shows  results  for  beams  all  of  length  to  diameter  ratio  approximately  11.2 
but  of  a  wide  range  of  hardnesses.  It  is  clear  that  the  strength  of  the  beams 
has  influenced  the  amount  of  bending  but  not  the  position  of  the  major  hinges. 

Strain  distributions  from  the  tension  side  of  all  beams  in  Fig.  1(a) 
and  one  from  Fig.  2  are  shown  in  Fig.  3  together  with  strain  profiles  from  the 
computer  simulations  to  be  discussed. 
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4.  DISCUSSION 


(a)  Rigid-Plastic  Approach 

The  development  and  motion  of  the  hinges  in  impacted  beams  is 
readily  described  using  the  rigid  plastic  approach  developed  by  Lee  and 
Symonds  [3]  for  central  impacts  on  a  free-free  beam.  In  Lee  and  Symonds'  case 
a  beam  of  length  21  is  in$>acted  in  the  centre,  at  which  position  a  plastic 
hinge  develops  if  the  load  is  sufficient  for  the  limit  moment  of  the  beam  to 
be  exceeded.  To  treat  the  end  impact  situation  we  use  the  same  geometry  but 
allow  the  centre  of  the  beam  to  have  a  limit  moment  equal  to  zero;  thus  the 
beam  yields  at  any  load  greater  than  zero  and  both  rotates  about  the  centre 
and  accelerates  in  the  direction  of  the  impulse.  It  is  then  possible  to  focus 
on  the  right  hand  section  of  length  1  as  a  simulation  of  our  end  impact 
problem,  the  moment  at  both  ends  of  one  beam  being  zero  and  the  first  stage  of 
the  end  impact  problem  corresponds  to  the  second  stage  of  the  central  impact 
problem  of  Lee  and  Symonds  [3]. 

The  situation  is  depicted  in  Fig.  4  where  the  impulsive  force  P  is 
applied  to  a  beam  of  length  21.  At  low  impulsive  forces  each  half  of  the  beam 
will  behave  in  a  rigid  manner  with  a  linear  acceleration  aQ  at  the  impact  site 
and  an  angular  acceleration  o  ,  being  positive  clockwise  for  the  right  hand 
half  which  is  to  be  considered.  Using  the  same  notation  and  method  as  Lee  and 
Symonds  [3]  aQ  and  aQ  are  given  as 


ao 

3 

2P/ml 

(2a) 

a 

o 

3 

3P/ml2 

(2b) 

where  m  is  the  mass  per  unit  length 

of 

the  beam. 

The  mean  acceleration,  a,  of 

a  length  of  beam  distance  x  from  the  impact  site  and  of  mass  (l-x)m  is  given 
by 


a 

a  -  a - 2.  (x  +  i) 

O  2 


(3) 


Hence  the  shear  force  (Vx)  and  bending  moment  (Mx)  at  x  are  given  by 


V  »  m(l  -  x)a  +  2  (x2  -  i2)  a  (4a) 

X  0  2  O 

ma  ma  _  la  .2 

MX  -  “T  x  -fX+"l  <*o  >  x  +-^?  (2V  -  ^o5  (4b) 
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Equating  Vx  to  zero  and  substituting  from  (2)  for  aQ  and  a  shows  that  the 
maximum  bending  moment  in  the  beam  occurs  where  x  *  1/3,  wRich  is  where 
initial  yielding  of  the  beam  can  be  expected  as  P  increased,  and  confirms  the 
result  of  Zener  and  Peterson  [1]. 

If  the  load  is  sufficient  to  exceed  yield  then  there  are  two 
segments  in  the  length  1  of  the  right  half  of  the  beam  joined  by  a  plastic 
hinge.  The  left  segment  has  a  linear  acceleration  aQ  at  the  inpact  site  and 
an  angular  acceleration  a  ,  and  the  right  hand  segment  has  a  linear 
acceleration  a^  at  the  en§  of  the  beam  and  an  angular  acceleration  a,. 
Corresponding  angular  velocities  in  the 
expressed  in  terms  of  the  limit  moment. 


M  “ 


The  position  of  the  hinge,  xh,  is  given 

5  = 


sections  are  and  The  load  is 

Mq,  of  the  beam  using  the  parameter^, 


PI 

M 


(5) 


dimensionlessly  as  £  by 


Vi 


(6) 


where  £  may  be  finite.  The  equations  of  motion  of  the  two  segments  become 
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ml' 

M 

o 
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3  3 

mlJ  (1  -  £)J 


M 


ml" 


V  _  6  6 

C  £2  (1  -  l)2 


(7a) 


(7b) 


(7c) 


(7d) 


-  £  (U)Q  -  u>,  ) 


(7e) 


These  equations  are  only  slight  modifications  of  the  equations  given  by  Lee 
and  Symonds  13]  for  the  central  impact  problem.  With  equations  (5),  (6)  and 
(7)  and  the  velocity  and  acceleration  formulas  established  by  Lee  and  Symonds 
[3]  it  is  possible  to  calculate  numerically  the  positions  of  hinges  and  the 
profile  of  a  beam  as  a  function  of  time  for  the  loading  function  P.  In 
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comparing  the  calculations  with  end  impact  experiments  it  is  noted  that  P  is 
obtained  by  doubling  the  load  applied  to  the  end  of  a  single  beam  in  the 
experiment. 


Numerical  evaluation  of  the  equations  shows  that  the  beam  rotates  as 
a  rigid  body  as  the  load  increases  until  a  hinge  develops  one  third  rod  length 
from  the  inpact  site.  With  further  increase  in  load  the  position  of  the  hinge 
moves  either  towards  the  impact  site  for  slender  beams  or  towards  the  centre 
for  beams  of  low  length  to  diameter  ratio.  If  the  load  P  is  sufficient  to 
cause  shear  yielding  of  the  end  of  the  beam  then  a  limit  is  reached  [6,7,8] 
because  shear  sliding  occurs  on  the  end  at  constant  stress.  In  this  case  the 
beam  is  effectively  subjected  to  a  constant  end  load  and  the  hinge  becomes 
stationary.  If  the  load  decreases  to  zero  the  plastic  hinge  moves  towards  the 
beam  centre,  the  beam  continuing  to  bend  as  the  hinge  moves  until  the  angular 
velocities  on  either  side  of  the  hinge  become  equal  at  the  time  the  hinge 
reaches  the  centre. 


For  a  circular  beam  the 
exceeded  is  given  by 


load,  P,  at  which  the  shear  yield  stress  is 


_P 

2 


TT 

777 


Y  D 


(8) 


The  limit  moment,  Mq  for  a  circular  beam  is 


M  =.  1/6  Y  D3 

o  ' 


(9) 


Substituting  (8)  and  (9)  in  (5)  and  putting  u>  =  u^and  evaluating  equation 
(7e)  gives  the  hinge  position  as  a  function  op  length  to  diameter  ratio  of  the 
beam,  as  shown  in  Fig  5.  At  times  greater  than  zero,  u>  >  necessarily  and 
if  the  load  exceeds  the  shear  yield  stress  P  will  continue  to  be  given  by 
equation  (8),  hence  y  is  constant  and  the  sign  of  £  can  be  evaluated  using 
equation^ 7e)  for  hinge  positions  to  either  side  of  the  initial  value.  The 
sign  of  £  is  such  as  to  return  the  hinge  to  its  initial  position,  hence 
whilst  y  remains  constant,  bending  will  continue  about  the  initial  hinge 
position.  The  bending  moment  at  yield  of  the  outer  fibres  is  3n/16  times  the 
limit  moment  given  by  equation  (9)  and  hinge  positions  evaluated  using  this 
yield  moment  are  also  given  in  Fig  5.  For  long  length  to  diameter  ratio  beams 
Fig  5  shows  the  hinge  positions  approach  limits  in  terms  of  beam  diameter 
whichever  moment  value  is  used.  The  value  approached  for  the  limit  moment 
case,  equation  (9),  is  approximately  1.1  beam  diameters,  and  this  agrees  with 
the  value  found  by  adapting  Symonds  [7]  later  work  on  beams  of  infinite 
length. 


The  major  hinge  near  the  impact  site  was  observed  to  be  in  an 
approximately  constant  position  independent  of  the  beam  strength  or  length  to 
diameter  ratio,  Figs.  1  and  2,  which  is  expected,  based  on  the  above  reasoning 
for  the  situation  in  which  the  shear  yield  stress  is  exceeded,  and  for  the 
length  to  diameter  ratios  used.  However,  the  hinge  position  in  all  cases  is 
between  two  and  four  beam  diameters  from  the  centre  of  impact:  the  average  is 
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very  close  to  three.  This  is  much  larger  than  the  1.1  diameters  expected 
using  the  limit  moment  in  the  theory  described  above  and  the  departure  is 
still  greater  if  the  yield  moment  is  considered,  Fig.  5.  Quasi-static  four 
point  bend  tests  on  some  of  the  beams  showed  that  for  these  beams  and  angles 
of  bend  the  yield  moment  was  a  better  estimate  of  the  bend  strength.  If  the 
inner  edge  of  the  impact  site  is  used  instead  of  the  centre  of  impact  then 
there  is  a  slight  improvement  but  still  the  major  hinge  is  too  far  from  the 
end  of  the  beam,  the  model  has  a  limitation  in  that  it  requires  bending  to 
occur  at  a  point  whereas  in  practice  one  would  expect  the  width  of  the  hinge 
to  be  at  least  of  the  order  of  the  beam  diameter,  which  is  the  order  of  the 
prediction  of  hinge  distance  from  the  inpact  site. 


Estimates  of  projectile  contact  time  from  computer  simulations 
(below)  were  between  20  and  30  microseconds.  Yield  stresses  for  equations  (8) 
and  (9)  were  estimated  as  2.9  times  the  Vickers  hardness  of  the  beams. 

Allowing  a  constant  load,  exceeding  the  shear  yield  stress,  r  the  end  of  the 
beam  for  these  contact  times  gave  bend  angles  which  were  toe  irje  but  the 
deflection  of  the  impact  point  was  predicted  to  well  within  e  correct  order 
of  magnitude  and  also  varied  in  the  correct  manner  with  len<  ,  to  diameter 
ratio  using  the  rigid  plastic  approach.  More  precise  ealeu  ions  would 
require  an  exact  knowledge  of  the  loading  history. 


(b)  Computer  Simulation 

Computer  simulation  offers  the  opportunity  to  study  the  detailed 
development  of  the  hinges  with  time  and  at  the  same  time  the  loading  can  be 
readily  handled.  The  computer  code  used  was  EPIC-2  which  is  a  finite  element 
code  in  two  dimensions.  Thus  the  simulation  is  in  plane  strain  although  the 
event  being  simulated  is  three  dimensional  and  exact  agreement  should  not  be 
anticipated.  Most  simulations  used  eight  triangular  elements  across  the  beam 
depth,  although  one  problem  was  also  run  with  twelve  elements  across  to 
determine  whether  grid  refinement  had  any  significant  effect.  Failure  strains 
for  the  beam  were  made  large  so  that  only  plastic  flow  was  allowed,  no 
fracture  having  been  observed  in  experiments. 


The  effective  strains  in  surface  elements  are  plotted  in  Fig.  3  for 
comparison  with  experimental  measurements.  The  magnitudes  of  the  strains 
agree  well  with  the  experimental  results  except  for  the  case  of  Fig.  3(a),  and 
as  noted  previously  in  reference  to  Fig.  1(a),  the  strains  in  this  beam  were 
unusually  low.  The  maxima  in  the  strain  distributions  from  the  simulation  are 
in  approximately  the  correct  position  and  in  changing  from  a  soft  to  a  hard 
beam  at  the  same  length  to  diameter  ratio  the  main  peak  is  observed  to  broaden 
and  its  magnitude  to  decrease  in  both  the  experiment  and  the  simulation.  Figs. 
3(c)  and  (e).  The  simulation  gives  two  distinct  peaks  in  strain,  one  near  the 
impact  site  and  one  near  the  centre  of  the  beam,  and  some  minor  peaks  between 
these.  The  mean  position  of  the  major  hinge  from  the  code  solutions  was  2.5 
diameters  from  the  impact  site  for  those  cases  simulated.  The  actual  position 
varied  slightly  with  the  length  to  diameter  ratio  and  strength  of  the  beam  but 
whether  or  not  this  is  a  true  variation  is  not  evident  from  the  solutions  as 
the  uncertainty  due  to  grid  size,  (+  one  quarter  beam  diameter),  is  of  the 
order  of  the  variation  in  position.  The  use  of  a  finer  grid  size  only 
slightly  modifies  the  observed  strain  distribution,  Fig.  3(c).  As  can  be 
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expected  from  the  description  of  the  strain  distributions,  the  final  profiles 
from  the  computer  simulation  were  very  close  to  those  observed  in  experiments. 


A  plot  of  a  typical  strain  distribution  across  a  beam  is  shown  in 
Fig.  6(a).  Adjacent  to  the  impact  site  the  plastic  strain  contours  go  almost 
straight  across  the  beam  which  confirms  the  estimate  that  the  loading  was 
sufficient  to  cause  shear  yielding  of  the  end  of  the  beam.  Further  along  the 
beam  the  symmetrical  nature  of  the  strain  profile  confirms  that  the  problem  is 
essentially  one  of  simple  bending  as  described  in  the  Lee  and  Symonds  [3] 
model.  The  development  of  the  strain  in  the  outer  fibres  as  a  function  of 
time  is  shown  in  Fig.  6(b).  The  rise  in  strain  towards  the  impact  site  is 
associated  with  localized  impact  damage.  Up  to  30  microseconds,  during  which 
time  the  load  is  applied,  a  peak  develops  around  one  point  in  the  beam 
although  the  effect  of  a  finite  width  is  evident  as  the  peak  is  broad.  After 
30  microseconds  the  projectile  had  left  the  beam  but  plastic  bending  continues 
to  occur  because  of  the  momentum  of  the  beam,  and  the  hinge  moves  towards  the 
beam  centre,  as  described  above  using  the  rigid  plastic  model.  When  the  hinge 
reaches  the  centre  bending  continues  and  develops  another  large  strain  peak, 
contrary  to  what  was  expected  from  the  rigid-plastic  model.  The  velocity 
profiles  at  30,  50,  70  and  100  microseconds  are  shown  in  Figs.  6(c)  to  (f)  for 
the  case  shown.  At  30  microseconds  the  situation  is  far  from  the  idealization 
of  two  rigid  sections  bending  about  a  hinge.  At  this  stage,  however  the  major 
hinge  is  the  only  place  at  which  significant  plastic  deformation  is 
occurring.  Whilst  the  load  has  been  removed  at  50  microseconds  the  velocity 
distribution  is  still  very  complex  and  at  70  microseconds  a  comparison  of 
velocity  vectors  at  equal  distances  either  side  of  the  point  of  rotation  shews 
that  bending  continues  to  occur.  At  100  microseconds  the  beam  is  essentially 
rotating. 


It  is  clear  from  these  velocity  pi  5files  that  although  most  plastic 
strain  is  associated  with  the  type  of  deformation  described  in  the  simple 
rigid-plastic  model,  both  during  loading  and  as  the  hinge  moves  to  the  centre 
after  load  removal,  in  the  experimental  situation  neither  segment  of  the  beam 
is  near  stress  equilibrium  at  any  stage.  If  the  impact  site  is  still  moving 
too  fast  and  the  distal  end  of  the  beam  too  slowly  when  the  hinge  reaches  the 
centre,  then  each  half  of  the  beam  can  be  elastically  bent  to  the  yield 
stress.  Locally,-  the  bending  may  cease  instantaneously  as  the  hinge  reaches 
the  centre  but  the  elastic  bending  in  the  two  halves  of  the  beam  will  then  add 
to  exceed  the  yield  stress  causing  further  strain  at  the  centre.  The  relief 
of  elastic  bending  of  the  segments  then  provides  a  possible  explanation  of  the 
strain  peak  near  the  centre  of  the  beams. 


5.  CONCLUSION 

Dynamic  impact  experiments  were  carried  out  on  the  end  of  free-free 
beams  of  circular  section  and  the  positions  of  hinges  and  the  overall  strain 
distributions  compared  with  rigid-plastic  model  and  computer  code 
predictions.  The  basic  deformation  mechanism  described  by  the  model  is 
confirmed,  but  it  predicts  a  hinge  position  too  close  to  the  impact  site 
because  it  does  not  take  account  of  the  finite  width  of  the  hinge.  In  the 
experiments  a  further  major  bend  was  observed  near  the  centre  of  the  beams, 
probably  associated  with  the  relief  of  elastic  stress  in  the  two  segments  of 
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the  bean,  a  factor  which  also  can  not  be  included  in  the  rigid  plastic 
model.  The  rigid-plastic  model  provides  a  clear  explanation  for  the 
development  of  a  stationary  hinge  whilst  the  beam  is  loaded  and  for  the 
movement  of  the  hinge  when  the  load  is  removed.  Hie  use  of  a  finite  element 
computer  code  allowed  the  loading  cycle  on  the  beam  to  be  simulated,  giving 
extremely  good  agreement  with  the  experimental  strain  distributions, 
confirming  the  nature  of  the  deformation  as  essentially  a  bending  process. 
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FIG.  2 


Beams  of  length  to  diameter  ratio  approximately  11.2  and  of 
hardnesses,  top  to  bottom,  575HV10,  480HV1 0,  450HV10,  395HV10  and 
355HV10. 
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FIG.  3  Strain  distributions  for  the  tension  side  of  beams  of  hardness 
395HV10  and  of  length  of  diameter  ratios  (a)  7.3,  (b)  11.6,  (c) 
15.3  and  (d)  23.8  from  Fig.  1(a),  and  (e )  of  length  to  diameter 
ratio  11.2  and  hardness  575HV10  from  Fig.  2.  In  each  case  the 
histogram  represents  experimental  measurements  from  the  bent 
beams  and  the  continuous  line  the  results  of  simulations  using 
the  EPIC-2  computer  program.  Computer  simulations  were  done 
with  two  grid  sizes  for  case  (c). 


PIG.  4  (a )  Impulsive  force  P  applied  to  a  beam  of  length  2  1  as  in  Lee 
and  Symonds  (3]  but  the  limit  moment,  Mp,  is  set  to  zero  at 
the  central  hinge  so  that  the  beam  bends  for  any  finite  load. 

(b)  Motion  of  rigid  beam  of  length  1  at  loads  below  that 
sufficient  to  cause  yielding.  The  load  P/2  applies  to  the 
beam  length  1  which  experiences  a  linear  acceleration  a0  at 
the  inpact  site  and  an  angular  acceleration  <*q. 

(c)  At  loads  sufficient  for  the  limit  movement  Mp  to  be 
exceeded,  two  rigid  segments  are  joined  by  a  plastic  hinge  the 
limit  moment  being  given  by  MQ  in  equation  (9).  'Hie  linear 
accelerations  of  the  inpact  and  distal  ends  are  aQ  and  a., 
respectively  and  the  angular  accelerations  of  the 
corresponding  segments  are  aQ  and  <a^  • 


Mo  *  Limit  Moment 
-  Mo  ■  Yield  Moment 


20  30 


FIG.  5  Hinge  position  divided  by  beam  diameter  (Xh/D)  as  a  function  of 
length  to  diameter  ratio  (L/D)  for  the  case  where  shear  yielding 
occurs  at  the  impact  site.  Two  situations  are  shown,  in  one  case 
the  beam  bends  plastically  at  the  limit  moment  (Equation  (9))  and 
in  the  other  it  bends  plastically  at  the  yield  moment. 
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Strain  distribution  from  the  computer  simulation  in  a  beam  of 
length  to  diameter  ratio  11 #6  and  of  the  same  strength  as  the 
beams  of  Pig.  1  (a).  The  strain  contours  have  values  of 
natural  strain  0.01,  0.02,  0,04,  0.06  and  0.07  and  the  strain 
increases  on  moving  out  from  the  neutral  axis  of  the  beam. 


Values  of  effective  strain  at  the  surface  as  a  function  of 
fractional  position  along  the  beam  for  the  same  problem  as 
above  at  the  times  indicated  on  the  curves,  in  microseconds. 
Only  that  part  of  the  curve  which  has  changed  with  increase  in 
time  is  shown. 


s  %  \ 


(e) 


(f )  Velocity  vectors  along  the  beam  for  the  case  shown  above 
at  times  (c)  30  microseconds,  (d)  50  microseconds,  (e)  70 
micro-seconds  and  (f )  100  microseconds. 
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